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应 用 修订 的 Shuttleworth-Wallace 模型 对 
半 干 旱 区 帮 腊 玉米 蒸 散 的 研究 
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$2 要 蒸 散发 (ET) 包 括 植物 蒸腾 (T) 和 土壤 蒸发 (E)， 在 维持 全 球 能 量 平 衙 和 气候 调节 中 起 关键 作用 。 量 化 蒸 
散发 及 其 组 分 在 准确 预报 气候 对 生态 系统 碳水 通 量 和 能 量 的 响应 中 至 关 重 要 。 基 于 兰州 大 学 半 干 旱 区 农业 生 
态 系统 试验 站 2014 年 玉米 生长 季 涡 度 相 关 仪 的 观测 结果 , 利用 修订 后 的 Shuttleworth-Wallace 模型 (S-W 模型 ) 
对 覆 膜 玉米 田 的 蒸 散 发 进行 模拟 ， 利 用 实测 值 对 模拟 结果 进行 验证 ， 对 蒸 散发 及 其 组 分 的 影响 因素 和 敏感 性 
进行 分 析 。 ERRA: S-W 模型 对 覆 膜 玉米 农田 蒸 散发 的 模拟 结果 在 日 蒸 散 量 大 于 2 mm':d 的 晴天 和 时 晴 时 
云天 气 较 好 ， 阴 雨天 气 模拟 结果 较 差 ， 且 模型 模拟 结果 较 涡 度 相 关 的 实测 值 偏 高 。E/ET 在 一 天 内 的 变化 为 单 
峰 曲 线 ， 在 生长 季 尺 度 上 ,在 玉米 快速 生长 期 呈 下 降 趋势 ,在 之 后 的 时 间 基 本 保持 不 变 。. 履 膜 玉 米 农田 的 E/ET 
在 日 时 间 尺 度 的 变化 主要 受气 孔 导 度 影响 ,在 生长 季 尺 度 主要 受 叶 面积 指数 和 土壤 含水 量 的 共同 影响 。 敏 感 
性 分 析 表 明 ， 蒸 散发 及 其 组 分 对 作物 冠 层 高 度 与 参考 高 度 间 的 空气 动力 阻力 (ra9) 和 冠 层 内 边界 层 阻 力 (ra9) 均 较 
人 敏感， 对 作物 冠 层 阻力 (rs9 人 敏感 性 适中 ， 对 地 面 与 冠 层 间 的 空气 动力 阻力 (ra 和 下 垫 面 裸 土 表面 阻力 (rs9) 不 人 敏 
感 ， 在 应 用 S-W 模型 模拟 覆 膜 玉米 农田 蒸 散 量 时 ， 要 特别 注意 阻力 参数 r，、ra* 和 rs 的 合理 确定 。 
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Estimation of evapotranspiration in maize fields with ground mulching with 
plastic film in semi-arid areas using revised Shuttleworth-Wallace model 
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Abstract  Evapotranspiration (ET), which is comprised of evaporation from soil surface (E) and transpiration from vegetation 
(T), plays an important role in maintaining global energy balance and regulating climate. Quantifying partitioning of ET is 
particularly important for accurate prediction of climate response to ecosystem carbon, water and energy budgets. Using eddy 
covariance measurements in maize fields for the growing season at the Experiment Station of Agro-ecosystem in Semiarid 
Area (ESASA) of Lanzhou University, we ran the revised Shuttleworth-Wallace model (S-W model) partitioned 
evapotranspiration in maize fields under plastic film mulch conditions into evaporation and transpiration, validated the 
performance of the model for different time scales and under different weather conditions with measured eddy covariance 
values, analyzed the driving factors, and determined parameter sensitivity of ET and its components. The results suggested that 
the simulated ET in the study area was in good agreement with the measurements in both sunny and cloudy days, but the model 
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performed badly in rainy days. In diurnal timescale, the modified model performed well when ET was larger than 2 mmrd'! in 
both sunny and mostly cloudy days, and the ratios of simulated values by S-W model to measured values were close to the 
1:1 line. But the model slightly overestimated ET in rainy days. Solar radiation and temperature were key environmental 
factors influencing ET in maize fields under plastic film mulch that led to seasonal variations. In general, T accounted for a 
small fraction of ET in maize fields under plastic film mulching. Diurnal variation in E/ET followed a single-peak curve, the 
low point was observed at night. At seasonal timescale, E/ET decreased from 18% to 8% at jointing stage, and kept at 8% level 
at tasseling stage and filling stage. We compared our result with others and found that both film mulch and environmental 
factors affected the value of E/ET. The E/ET was controlled by canopy stomatal conductance at diurnal timescale, while at 
seasonal timescale, it was mainly controlled by leaf area index (LAI) and soil moisture content (0) which regulated transpired 
water from leaf stomata and evaporated water from bare soils. The sensitivity analysis showed that ET and its components 
were most sensitive to aerodynamic resistance from canopy to reference height (ra and bulk resistance of boundary layer (r^), 
and moderately sensitive to bulk resistance of canopy stomatal (r.^), and insensitive to aerodynamic resistance from soil to 
canopy (rą) and soil surface resistance (rē). It was suggested that ET and its components were more sensitive to parameters 
related to canopy. Therefore, it was much more important to determine resistance parameters of ra’, r,^ and rs when simulating 


ET in maize fields under plastic film mulch using S-W model. 
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TREHAKBIUR Gu FRTETIARS TU dD3RIAUR, ÆtR- 
植被 -大 气 连 续 体 (SPAC) 水 量 转换 和 能 量 转 移 的 重 
要 组 成 部 分 品 ， 多 年 来 一 直 是 气象 学 、 植 物 学 、 水 文 
学 、 生 态 学 、 地 理学 、 土 壤 学 等 众多 学 科 关 注 的 焦 
点 问题 之 一 上 ]。 作 物 蒸 散 与 气象 因子 、 土 壤 水 分 状 
况 、 作 物种 类 及 其 生理 特性 等 多 种 因素 相互 影响 ， 
在 一 定 程度 上 制约 着 小 尺度 微 气象 状况 和 大 范围 气 
候 变 化 ”在 中 国 西北 干旱 半 干 旱 区 ,气候 干燥 、 
降水 稀少 、 共 发 作用 强烈 ， 使 得 水 资源 成 为 该 区 域 
AERA MOO, 因此, 了 解 西北 地 区 作物 燕 散 发 
及 水 分 运 移 规律 ,不 仅 可 以 为 该 区 域 能 量 平衡 和 水 
分 循环 的 研究 提供 理论 依据 ， 还 对 确定 该 区 域 水 分 
利用 效率 和 水 资源 分 配 起 着 关键 作用 。 

在 众多 计算 蒸 散发 的 模型 中 ，Penman-Menteith 
(P-M) 模 型 无 疑 是 一 个 基准 公式 ， 该 模型 在 1998 年 
被 联合 国 粮 食 及 农业 组 织 (FAO) 推 荐 为 计算 参考 作 
物 蒸 散 量 的 标准 方法 02。1985 年 Shuttleworth 和 
Wallace! P-M 公式 的 基础 上 ， 将 研究 对 象 细 分 为 
作物 冠 层 和 冠 层 下 地 表 两 部 分 ， 建 立 了 主要 用 于 稀 
EE RB) Shuttleworth-Wallace(S-W) 双 源 蒸 
散 模 型 。 与 P-M 模型 相 比 ，S-W 模型 将 作物 蒸腾 和 
土壤 蒸发 分 开 计 算 , 考虑 了 作物 与 裸 土 的 阻力 差异 ， 
提高 了 蒸 散 发 的 计算 精度 。 已 有 学 者 将 两 种 模型 的 
模拟 值 进行 了 上 比较 , 如 Stannard HE 3E F E Ac 83 
自然 植被 状态 下 对 P-M 模型 、S-W 模型 和 Priestley- 
Taylor 修正 模型 进行 了 对 比 , ARRA, ENE 
被 状况 下 S-W 模型 的 模拟 结果 优 于 P-M 模型 ; Kato 
SEIDIEESS P-M 和 S-W 模型 对 稀 玖 植被 的 适用 性 ， 
结果 表明 S-W 模型 的 均 方 根 误差 远 远 小 于 P-M 模型 ; 
Fisher 等 9 在 对 两 个 模型 的 比较 中 也 得 出 了 相同 结 
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论 。 由 于 S-W 模型 中 参数 众多 ,尤其 是 几 个 阻力 参 
数 的 取 值 对 模型 影响 很 大 。Leuning 77, Sellers 
等 05 和 Zhu 等 2 分 别 引 入 饱和 水 气压 差 和 土壤 含水 
量 对 气孔 导 度 进行 优化 ; Anadranistakis EPORA E 
壤 含 水 量 优化 了 土壤 表面 阻力 参数 ; Zhu EPN 
叶 斯 方法 引入 到 阻力 参数 的 计算 中 , 提高 了 模型 的 
准确 性 ; Li 等 中 考虑 了 地 膜 覆盖 对 蒸 散 发 的 影响 ， 
对 S-W 模型 进行 了 优化 ， 结 果 发 现 优 化 后 模型 的 模 
拟 结果 优 于 传统 S-W 模型 ， 计 算 的 土壤 蒸发 更 接近 
于 实测 值 。 

玉米 (Zea mays L.) 是 世界 上 分 布 最 广泛 的 粮食 作 
物 之 一 ,种 植 面积 仅 次 于 小 麦 (TYiticum aestivum. L.) 
和 水 稳 (Oryza sativa L.) 而 居 第 3 位 ， 是 我 国 西北 干旱 
半 干 旱地 区 的 主要 农作物 ,也 是 当地 农民 的 主要 经 
济 来 源 。 地 膜 覆 盖 作 为 一 项 增产 措施 , 可 以 有 效 减少 
土壤 蒸发 二 近 些 年 在 我 国 北方 地 区 农业 生产 中 得 
到 了 广泛 应 用 呈 。 以 往 关 于 农田 生态 系统 蒸 散 发 的 研 
究 鲜 有 将 地 膜 覆 盖 的 影响 考虑 在 内 。 将 修订 的 S-W 模 
型 应 用 到 黄土 高 原 半 干 旱 区 覆 膜 玉米 农田 茹 散发 的 
研究 不 仅 可 以 检验 该 模型 的 适用 性 , 还 是 制定 合理 
的 栽培 和 管理 措施 的 基础 ， 对 提高 水 分 利用 效率 、 增 
加 当地 经 济 有 重要 意义 。 本 文 在 2014 年 生长 季 观 测 
资料 的 基础 上 , 用 修订 的 S-W 模型 对 覆 膜 玉米 农田 
蒸 散 发 进行 模拟 并 验证 ， 旨 在 全 面 研究 蒸 散发 规律 并 
探究 蒸 散 发 及 其 组 分 的 影响 因素 和 参数 敏感 性 ,为 确 
定 履 膜 玉米 合理 的 灌溉 制度 提供 一 定 的 理论 指导 。 


1 研究 区 概况 与 研究 方法 


11 ”研究 区 概况 
试验 田 位 于 兰州 大 学 半 干 旱 区 农业 生态 系统 试 
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验 站 (35.42°N,105.24°E, HHR 1951 m), 研究 区 属于 
黄土 高 原 丘 陵 沟 至 区 ,气候 类 型 为 典型 的 暖 温带 半 
干旱 季风 气候 ,年 平均 气温 7.9 C, 年 平均 降水 量 
340 mm, 主要 集中 在 7 一 9 月 , 年 潜在 蒸发 量 1 800 mm, 
无 霜 期 约 136~186 d, 日 夜 温差 大 ,光照 充足 。 冬 季 
干旱 多 风 ， 夏季 炎热 ， 地 表 水 稀少 ,地 下 水 位 深 。 区 
内 主要 包括 内 山 旱 地 、 拓 旱 地 以 及 川 旱地 ， 耕 地 严 
重 缺 水 ,粮食 作 物 以 玉米 为 主 ， 主 要 气象 灾害 是 干 
旱 。 试 验 田 为 常规 种 植 农田 ， 玉 米 行距 和 株距 分 别 
7J 0.5 m 和 0.4 m 本 研究 所 用 玉米 在 4 月 底 播种 (种 
植 时 铺设 地 膜 )，10 月 初 收获 , 7. 8 月 玉米 处 于 拔节 
期 、 抽 雄 期 和 灌浆 期 ， 此 时 玉米 蒸腾 作用 强烈 ， 在 此 
期 间 玉 米 农 田 保持 在 自然 状态 , TATEM. EX 
种 植 时 覆 膜 面积 占 农田 面积 的 65%， 随 着 玉米 的 出 
苗 、 生 长 和 天 气 等 外 部 因素 的 影响 ， 覆 膜 面 积 的 比 
例 降 低 至 60% 左 右 。 
1.2 ”数据 来 源 
1.2.1 ”农田 茸 散 数据 观测 

农田 蒸 散 发 用 涡 度 相关 仪 观测 ， 观 测 期 为 2014 
年 6 月 29 日 一 9 月 10 H, 在 长 势 较 好 的 玉米 农田 中 
央 架 设 通 量 观测 塔 ， 涡 度 相 关 仪 安装 在 通 量 观测 塔 
上 3 nm 高 度 , 包括 三 维 超声 风速 仪 (CAST3，Campbell， 
USA) 和 CO;/H;O 红外 气体 分 析 仪 (Li-7500，LiCor 
Inc.，USA)。 取 样 频率 为 10 Hz， 通 过 数据 采集 器 
(CR5000, Campbell, USA) 采 集 并 按 30 min 计算 平均 
值 进 行 存 储 ， 数 据 包括 感 热 通 量 、 潜 热 通 量 、CO， 
和 水 汽 通 量 等 。 另 外 , 净 辐 射 仪 CNRI(Kipp&Zonen, 
Netherlands) 也 安装 在 塔 上 3 m AE, 在 净 辐 射 仪 正 
下 方 土壤 5 cm 深 处 , 分 4 个 方位 ( 东 、 南 、 西 、 北 ) 埋 
设 热 通 量 板 , 土壤 热 通 量 依 据 通 量 板 数 据 平均 求 得 。 

使 用 涡 度 相关 数据 时 需 对 数据 进行 校正 处 理 ， 
本 研究 使 用 Eddypro 5.3.0 软件 对 数据 进行 处 理 ， 对 
处 理 过 的 数据 依据 以 下 原则 进行 剔除 : 1) 降 水 以 及 
前 后 1 h 数据 ; 2) 传 感 器 状态 异常 数据 ; 3) 明 显 超出 
物理 意义 的 数据 。 对 缺失 和 剔除 数据 ,通过 对 灌流 
通 量 日 变化 绘图 分 析 并 用 内 插 法 进行 插 补 中。 利用 
FSAM 模型 计算 出 该 区 域 在 玉米 生长 季 90% 的 通 量 
来 源 于 迎风 向 0~150 m， 垂 直 于 迎风 向 -45~45 m’, 
在 此 区 域 均 种 植 玉米 , 利用 EBR 方法 计算 出 该 区 域 
7 一 8 月 的 能 量 闭合 度 为 0.905, 说 明 该 数据 质量 较 
好 。 使 用 处 理 过 的 潜 热 数据 ， 将 其 转化 为 农田 蒸 散 
发 (ET) 数 据 。 
1.2.2 土壤 水 分 及 气象 数据 观测 

埋设 ECH2O 传感器 探头 于 样 地 地 表 以 下 10 cm 
深 处 ， 对 土壤 体积 含水 量 进行 连续 监测 ， 监 测 日 期 


为 2014 ££ 6 H 12 H—9 H 10 A, 设置 数据 记录 间 
隔 时 间 为 0.5 php。 常规 气象 数据 采用 试验 站 内 通 量 观 
测 塔 附近 的 Watchdog 2700 小 型 自动 气象 站 测量 结 
果 , 包括 太阳 辐射 、 大 气压 、 温 度 、 湿 度 、 风 向 、 
风速 和 降水 等 。 根 据 气 温和 相对 湿度 计算 饱和 水 汽 
压 亏 缺 。 
1.2.3” 叶 面积 指数 的 获取 
在 观测 期 每 隔 1 周 左右 在 试验 田 随 机 选取 8 株 
玉米 ,测定 玉米 株 高 、 每 株 玉米 叶片 数 及 各 叶片 的 
长 、 宽 , 用 下 面 公式 计算 叶 面 积 指数 (LAD): 
a (1) 
RA: Li 为 每 株 玉米 各 叶片 的 叶 长 , Bij 为 最 大 叶 宽 ， 
n 为 第 j 株 的 总 叶片 数 , m 为 测定 株数 ，pa 为 种 植 密 
度 。 另 外 ,在 测量 期 间 每 隔 1 周 观测 1 次 气孔 导 度 ， 
挑选 典型 晴天 ， 从 7:00 一 19:00, 每 h FB Li-6400 测 
量 1 株 玉米 上 上、 中、 下 3 处 叶片 的 正 、 反 面 气孔 导 
度 和 光合 强度 等 数据 。 
1.3 ”修订 的 S-W 模型 介绍 
Li 等 所 修订 后 的 S-W 模 型 如 图 1 所 示 。 总 的 潜 热 
通 量 包括 来 自 冠 层 (4 刀 、 裸 土 4EJ 和 地 膜 下 土壤 的 洪 
热 4E,)。 覆 膜 状况 下 作物 的 潜 热 通 量 用 下 式 表 达 : 
AET = AT +(1- f, )AE, + f, AE, 
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S 


-1 
R +R R R R 
C, =——=——— -1+ —+(1- fna) + 2 6 
SUR | 及 ( Ing Aen (6) 
-1 
R, +R R R R 
C, 2—— dl 2 十 a 7 
jene | E O-R nal () 
E 
R +R R R R 
Cn == .;1+—+(1 34 a 8 
"TUR | R. (1— fm) R. fm zl (8) 
RS(A*y) +y (9) 
R =(A+7) R +yr (10) 
Rna -(A*y)r tyre (11) 
R, 2(^* yy? (12) 


式 中 : PM.、PMs 和 PMm 分 别 为 彭 曼 公 式 计 算 的 植被 


http://www.ecoagri.ac.cn 


第 5 期 


= EF: 应 用 修订 的 Shuttleworth-Wallace 模型 对 半 干 旱 区 帮 膜 玉米 蒸 散 的 研究 13X yA HETT] 


蒸腾 、 裸 土 蒸发 和 覆 膜 土壤 蒸发 (wm ，C.、Cs 和 
Can 分别 为 冠 层 阻力 系数 、 土 壤 表 面 和 地 膜 阻 力 系数 
(无 量 纲 )，/ 为 水 蒸发 的 汽化 潜 热 (MJ.kg P), 4 为 饱和 
水 汽 压 -温度 曲线 斜率 (kPaK-)，p 为 空气 密度 
(kg:m 2), es 和 e, 分 别 为 参考 高 度 (3 m) 的 饱和 水 汽 压 
和 实际 水 汽 压 (kPa), y 为 湿度 计 常 数 , 刀 为 覆 膜 面积 
与 农田 面积 的 比例 ,x 为 作物 冠 层 高 度 与 参考 高 度 
ig B3 Ze 377 [R73 (s:m ), rs 为 冠 层 内 边界 层 阻 力 
(sm 0)，rs 为 地 面 与 冠 层 间 的 空气 动力 阻力 (sm-)， 
x 为 作物 冠 层 阻力 (s-m!), rs 和 x" 分 别 为 下 垫 面 宰 
土 和 覆 膜 土壤 表面 阻力 (s.-m')。4 为 总 有 效能 量 , Ac. 
4 和 4nm 分 别 为 冠 层 、 裸 土 和 地 膜 接 受 的 有 效能 量 
(W-m?), 分 别 通过 下 式 计算 : 


A= R, - (1- fu)G, - fun (13) 

A, = R, — R, exp(-cLAI) (14) 

A, = Ras — G; = Ra exp(-cLAD) - G, (15) 
A, = Ras — Gn = R, exp(-cLAI) - G,, (16) 


式 中 : GACA AREMA FERRNA (W m’), 
KA 为 光 豪 减 消光 系数 , LAI 为 叶 面 积 指 数 (m”.m)。 


e 边界 Boundary Awy 阻抗 Resistance 


——e 传输 方向 Transmission direction 


参考 高 度 


us n AE | 


E 
7a 


Je ES E 


1 地 膜 覆 盖 玉 米田 修订 后 的 Shuttleworth-Wallace 
模型 示意 
Fig. 1 Schematic diagram of the revised Shuttleworth- 
Wallace model of plastic film mulched maize field 


在 方程 (2) 和 (3) 中 ， 参 数 1 和 4 通过 气象 数据 计算 
得 到 ,根据 Allen 等 "的 方法 计算 : 


4 = 2.501-0.002 361T, (17) 
1727T. 
4 098[0.610 8 a 
l Exp G 12373! 


A- 


: (18) 
(T, 4237.3) 
RELIE HRS, GGG, HH fà a ARE 
IEKA ERRE, BDAE,-0. MAA: 
AET sw — AT +(1- f, )AE, = C,PM, +(1- f, )C,PM, (19) 


-1 

R +R R R 
C =a e. 1+ —+(1 3 20 
, R. | R. ( 1I (20) 

-1 

R +R R R 
C, = — .41+ 一 二 (1 à 21 
, 及 | R. ( 1I (21) 
A-R,-G (22) 
A, = R, exp(-cLAI) - G (23) 


14 ”模型 阻力 参数 计算 

本 试验 针对 不 同 的 阻力 参数 选取 了 不 同 的 计算 
方法 。 

土壤 表面 阻力 系数 (x;) 根 据 Anadranistakis 等 中 
提出 的 计算 公式 进行 计算 : 


rŠ x nminf (0) ES T min (a ^ = 3 (24) 


式 中 : 0:23 EH Bl EZK (cm em); rsoin 是 土壤 水 分 
为 Gt 时 的 土壤 表面 阻力 ， 取 100 sm, 0 为 表层 土 
壤 合 水 量 (cm cem; 并 根据 Zhang 等 2 的 研究 确定 
ai 分别 为 2.5 和 1.5。 
作物 冠 层 阻力 参数 (9 可 以 用 气孔 导 度 (g\9) 的 
倒数 来 表示 。 gs* 的 计算 参照 Zhu 等 ("的 计算 方法 计 
算 公式 如 下 : 
ge = Emax n| Oh + Oso | 1 
Q, exp(-KgLAD) + Oso || 1+ D, /Dso 


vo 


(25) 
1 (0>0,) 
(0,«0«0,) (26) 
“0” (0-8) 

式 中 : gmax 为 冠 层 高 度 的 叶片 最 大 气孔 导 度 (ms 1); 
Ko 为 短波 辐射 的 消光 系数 ; 2 为 冠 层 顶部 可 见 光 辐 
射 的 通 量 密度 (W.m“); 2so 为 气孔 导 度 为 最 大 值 一 半 
时 的 可 见 光 辐 射 (W.m-); 万 ,为 参考 高 度 的 水 汽 压 亏 
缺 (kPa); D;o 为 气孔 导 度 为 最 大 值 一 半 时 的 水 汽 压 亏 
缺 (kPa); 94 为 凋 萎 时 的 土壤 舍 水 量 ，0, 为 作物 胁迫 开 
始 时 的 临界 含水 量 , 0, 取 饱和 合 水 量 0, 的 3/4。 

边界 层 阻力 参数 (x,*) 和 空气 动力 学 参数 (xs、 n), 
fk HE Shuttleworth 4l Gurney’! 在 1990 年 提出 的 公式 
进行 计算 。 
2 结果 与 分 析 
2.1 观测 期 玉米 生长 环境 分 析 

玉米 生长 期 间 环 境 因 子 的 变化 对 蒸 散发 及 其 组 
分 的 影响 很 大 。 观 测 的 2014 年 6 月 29 日 一 9 月 10 日 净 
辐射 、 气 温 、 水 汽 压 、 风 速 、 叶 面积 指数 和 土壤 合 
水 量 的 变化 如 图 2 所 示 。 从 6 月 29 日 开始 玉米 进入 快 
速生 长 期 ， 株 高 和 叶 面 积 指数 增加 较 快 ， 到 7 月 20 日 


/(0)- 
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VARA 


之 后 ， 叶 宽 和 株 高 变化 相对 变 缓 ，8 月 12 日 开始 玉米 
进入 灌浆 期 ， 玉 米 株 高 不 再 长 高 , 叶 面 积 指数 也 逐 
渐 减 小 。 在 观测 期 间 日 平均 气温 在 13~22 人间 变化 ， 
8 月 2 日 之 前 气温 变化 较 大 ， 变 化 幅度 超过 9 C, 而 8 
月 2 日 之 后 气温 较 低 且 变化 幅度 较 小 , 为 13~18 Co 

观测 期 间 气 温 整体 上 呈现 出 先 增加 后 减 小 的 趋势 。 净 辐 
射 的 变化 趋势 与 气温 相似 , 在 0.3~18 MI-m 7d ' [B S (E, 


25 


20 


ul 
Temperature (°C) 


N 


风速 
Wind speed (m:s ) 


Uu 


WERE 
Leaf area index (m^-m ?) 


N 


1 


06-29 07-19 08-08 08-28 


饱和 水 汽 压 差 
Vapor pressure deficit (kPa) 


净 辐 射 


土壤 含水 


在 玉米 快速 生长 的 7 月 份 , 净 辐 射 维持 在 较 高 水 平 。 
饱和 水 汽 压 差 在 7 月 份 波动 较 大 ， 最 大 值 为 1.35 kPa, 
从 8 月 初 开 始 减 小 。 在 观测 期 风速 呈现 出 先 增 加 后 
减 小 的 趋势 ， 日 平均 风速 为 0.5~2.5 mso 10 cm 深 
土壤 体积 含水 量变 化 幅度 较 大 ， 受 降雨 量 影响 ， 
在 0.17 cmi.cm“ 至 饱和 含水 量 之 间 变 化 ,饱和 含水 
量 确定 为 0.3 cm .cm ^, 
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日 期 (月 -日 ) Date (month-day) 
2 各 环境 因子 在 玉米 生长 季 内 的 变化 


Fig.2 Variations of environmental factors in growing season of maize 


2.0. S-W 模型 模拟 蒸 散发 与 实测 值 对 比 

为 了 研究 不 同 天 气 和 时 间 尺 度 的 玉米 蒸 散 ， 本 
文选 取 晴 天 、 时 晴 时 云 和 阴雨 天 3 种 典型 天 气 在 半 小 
时 尺度 上 对 玉米 蒸 散 进 行 了 验证 ,结果 如 图 3a、b、c 
所 示 。3 种 天 气 的 划分 依据 初 小 静 等 上 的 方法 。 从 图 
中 看 出 晴天 玉米 蒸 散 速率 最 高 可 达 0.5 mmh, A 
蒸 散 量 平均 在 3.53 mm.d-1，S-W 模 型 模拟 结果 与 涡 度 
相关 仪 实测 值 具有 较 好 的 相关 性 ， 二 者 拟 合 直线 的 
决定 系数 R* 为 0.917; 在 时 晴 时 云天 气 , 在 太阳 辐射 
较 强 时 ， 蒸 散 速 率 最 大 值 与 晴天 相近 ， 而 在 多 云天 
气 下 则 较 小 ， 为 0.2~0.3 mmh, ARAE RJ 


1.5~2.5 mm.d"!， 此 时 S$-W 模 型 模拟 值 与 实测 值 的 相 
关 性 也 较 好 ，R’ 为 0.869; EARRAS, EKA BOR 
率 明 显 降 低 ， 最 大 蒸 散 速率 平均 约 为 0.2 mmh”, E 
米 日 蒸 散 量 也 较 小 ,一般 为 0.3~1.0 mm-d^', 而 S-W 
模型 模拟 值 与 蒸 散 实测 值 的 拟 合 结果 相对 较 差 ,及 
为 0.733。 对 上 比 不 同时 间 的 模拟 结果 发 现在 晴天 和 时 
晴 时 云天 气 ，S-W 模 型 的 模拟 效果 更 好 ， 而 在 太阳 辐 
射 较 小 的 阴雨 天 气 ， 该 模型 模拟 结果 较 差 。 

日 时 间 尺 度 上 的 S-W 模 型 模拟 的 玉米 农田 蒸 散 
与 实测 值 的 比较 如 图 3d 所 示 。 相 关 性 分 析 表 明 , S-W 
模型 模拟 的 燕 散 与 涡 度 相关 仪 观测 的 燕 散 相关 系数 
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7- H679 


0.6 


a ”晴天 Sunny days 


模拟 值 Simulated value (mm) 


0 0.1 0.2 


0.6 


b ”时 晴 时 云 Mostly cloudy days 


0.4 


实测 值 Measured value (mm) 


3 不 同时 间 尺 度 和 天 和气 下 覆 膜 玉米 田 蒸 散发 ET) 模 拟 值 和 实际 值 的 比较 
Fig. 3 Comparison of simulated and measured evapotranspiration (ET) of plastic film-mulched maize field on different time 
scales and weather conditions 
a. b. c 为 半 小 时 尺度 , d 为 日 尺度 。Figure a, b and c are on the half-hour time scale, figure d is on the daily scale. 


为 0.909。 结 合 图 中 S-W 模 型 计算 结果 与 实测 值 的 关 
系 ， 可 以 发 现 该 模型 的 日 蒸 散 量 模拟 值 与 实测 值 的 
散 点 大 多 分 布 在 1 : 1 直线 以 上 ,模拟 结果 较 涡 度 相 
关 法 的 实测 值 偏 高 。 另 外 ,还 可 以 发 现 ,在 日 蒸 散 量 
大 于 2 mmd 时， 模拟 值 与 实测 值 的 点 较为 集中 地 
分 布 在 1 : 1 直线 附近 , 拟 合 效 果 较 好 ; mute HEX 
量 小 于 2 mm.d 的 范围 内 , S-W 模 型 的 模拟 值 与 实测 
值 的 点 偏离 1 : 1 直线 较 多 ， 拟 合 结果 较 差 。 
2.5 S-W 模型 模拟 土壤 蒸发 占 蒸 散发 比例 E/ET) 在 

不 同时 间 尺 度 上 的 变化 

利用 S-W 模 型 模拟 的 覆 膜 玉米 农田 土壤 蒸发 占 
共 散 发 比例 (E/ET) 日 变化 和 生长 季 内 变化 如 图 4a 和 
图 4b 所 示 。 从 图 4a 可 以 看 出 E/ET 在 一 天 内 的 变化 为 


E/ET 


0 
0:00 4:00 8:00 12:00 16:00 20:00 0:00 


时 刻 Time 


单 峰 曲线 , 与 净 辐 射 的 变化 相似 ，E/ET 在 夜晚 保持 
在 较 小 值 并 处 于 波动 之 中 ,因为 夜晚 土壤 蒸发 较 弱 ， 
峰值 则 出 现在 14:00 点 。 从 图 4b 可 知 在 拔节 期 随 着 玉 
米 的 快速 生长 E/ET 呈 快速 下 降 趋 势 ， 从 18% 降 至 8% 
左右 ; 在 抽雄 期 和 灌浆 期 ， 玉 米 农田 的 郁 闭 度 较 高 ， 
E/ET 变 化 较 小 ,在 8% 上 下 浮动 。Ferretti 等 652 对 科 罗 
拉 多 州 东北 部 矮 草 草原 的 蒸 散 发 进行 研究 ， 发 现 
E/ET 的 平均 值 为 7%; Sauer ERA E [Giycine 
max (L.) Merrill] 作 物 的 研究 中 发 现 当 植被 郁 闭 度 较 
高 (LAI>5 m^m 时 E/ET 小 于 8%:; Zhu 等 2 利用 S-W 
模型 对 甘肃 省 黑河 流域 的 玉米 农田 蒸 散 发 进行 了 研 
究 ， 模 拟 得 到 的 E/ET 为 13%。 本 研究 得 到 的 E/ET 比 
Zhu 的 结果 低 ， 这 可 能 是 因为 本 试验 的 玉米 农田 覆盖 


0.04 - 
06-30 07-14 07-28 08-11 08-25 09-08 
日 期 (月 -日 ) Date (month-day) 


4 ” 覆 腊 玉米 农田 土壤 蒸发 占 蒸 散 发 比例 (E/ET) 在 日 尺度 (a) 及 生长 季 尺 度 (b) 的 变化 
Fig.4 Variations of ratio of evaporation to evapotranspiration (E/ET) in plastic film-mulched maize field on diurnal (a) and 
seasonal (b) timescales 
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BA 


地 膜 ， 有 效 减少 了 裸 地 土壤 水 分 的 蒸发 进而 影响 E/ET。 
而 Ferretti 等 BA 和 Sauer 等 所 的 结果 与 本 研究 结果 相 
似 , 分 析 其 原因 ,除了 地 膜 外 ,不 同 植物 的 郁 闭 度 
也 可 能 对 E/ET 有 较 大 影响 。 为 了 研究 环境 因素 对 
E/ET 的 影响 , 本 试验 分 析 了 生长 季 内 叶 面 积 指数 、 


表 1 


土壤 合 水 量 和 气孔 导 度 与 B/ET 的 相关 关系 , 如 表 1 
所 示 。BE/ET 与 叶 面积 指数 和 和 气孔 导 度 呈 负 相关 ,与 
土壤 含水 量 呈 正 相 关 。 叶 面积 指数 和 气孔 导 度 主要 
是 通过 增加 叶片 蒸腾 所 占 的 比重 进而 影响 E/ET, E 
壤 合 水 量 则 通过 增加 土壤 蒸发 的 水 量 来 影响 E/ET。 


履 膜 玉 米 农田 土壤 蒸发 占 蒸 散发 比例 E/ET) 与 作物 和 环境 因素 的 拟 合 模型 


Table 1 Different models fitted for the ratio of evaporation to evapotranspiration (E/ET) and meteorological and crop factors 
of plastic film-mulched maize field 
因素 Factor 回归 方程 Regression equation R RSME n 
叶 面 积 指 数 Leaf area index E/ET=-0.38xLAI+0.201 0.581 0.016 65 
土壤 合 水 量 Soil moisture content E/ET-0.77x0—0.067 0.320 0.020 65 
气孔 导 度 Stomatal conductance E/ET=-0.989xgs+0.108 0.042 0.024 65 


2.4 ”模型 参数 敏感 性 分 析 

为 了 探究 蒸 散发 及 其 组 分 对 5 个 阻力 参数 的 敏 
感性 ， 本 试验 通过 计算 各 阻力 参数 变化 +10% 时 蒸 散 
发 及 其 组 分 的 变化 来 对 其 分 析 , 结果 如 表 2 所 示 。 从 
表 2 可 以 看 出 对 蒸 散 发 (ET) 影 响 最 大 的 参数 是 x,*， 当 
其 改变 10% 时 蒸 散发 改变 3.30%~3.55%; 其 次 是 re 
和 x“， 这 两 个 参数 增加 和 减 小 10% 时 蒸 散 发 分 别 减 
小 1.34%~1.74% 和 增加 1.40%~2.35%; Tür, Rer, x42 
散发 的 影响 最 小 ， 当 其 改变 10% 时 蒸 散发 改变 0.7%4 
左右 。 玉 米 蒸 腾 (T) 对 r* 和 六 "敏感 ， 当 这 两 个 参数 改 
变 10% 时 蒸腾 改变 1.32%~1.60%; 其 次 对 rs: 敏感 性 
适中 ; 对 和 ru 不 敏感 ， 当 rs* 和 7 变化 10% 时 蒸腾 改 


表 2 


变 不 到 0.2%。 蒙 发 (E) 对 re 和 ms 都 很 敏感 ， 当 me 和 7 
变化 10% 时 ,蒸发 改变 4.65%~5.39%; RE Ar Mrs 
也 较 敏 感 ， 当 其 改变 10% 时 蒸发 分 别 改变 2% 和 3% 左 
右 ; 而 当 " 改 变 10% 蒸 发 只 改变 0.31%~0.36%。 总 的 
KH, RRE, EKREM ERRERA r Mr 
感 ， 对 x“ 的 敏感 性 适中 ， 对 x 都 不 敏感 ， 这 说 明 冠 
层 的 阻力 系数 对 蒸 散 发 及 其 组 分 的 影响 更 大 。 李 俊 
等 2 在 用 S-W 模 型 和 P-M 模 型 对 麦田 蒸 散 发 的 敏感 
性 分 析 中 得 出 用 两 个 模型 估算 的 ET 均 对 冠 层 阻力 
最 敏感 ， 这 和 与 本 研究 结果 相 一 致 。 因 此 在 应 用 S-W 
模型 计算 农田 蒸 散 发 时 ， 应 注意 对 r*、re 和 六 的 合 
理 确 定 。 


履 膜 玉米 农田 燕 散发 及 其 组 分 对 阻力 参数 变化 +10% 的 敏感 性 


Table2 Sensitivity of simulated evapotranspiration and its components of plastic film-mulched maize field to the +10% changes in 
resistance parameters 


阻力 参数 蒸腾 Transpiration (%) 蒸发 Evaporation (%) 蒸 散 发 Evapotranspiration (96) 
Resistance parameter +10% -10% +10% -10% +10% -10% 
r -1.32 1.38 4.98 —-5.14 3.30 -3.55 
ry 0.12 —0.14 —0.31 0.36 —-0.71 0.76 
ra 一 1.44 1.60 —4.65 5.39 一 1.74 2.35 
ry 一 1.00 1.04 2.19 -2.25 一 1.34 1.40 
rs 0.07 —0.07 —3.08 3.33 0.71 —0.82 


3 Wit 
3.1 “影响 覆 膜 农田 蒸 散 发 的 因素 

地 膜 的 主要 作用 是 垃 温 、 保 增 ， 将 无 效 的 土壤 茹 
发 转化 为 有 效 的 作物 蒸腾 。 由 于 地 膜 的 阻隔 作用 ， 除 
了 少量 水 分 通过 地 膜 开 孔 蒸 发 到 空气 中 ， 其 余 水 分 只 
能 通过 作物 根部 吸水 而 释放 到 大 气 中 , 这样 就 增加 了 
作物 的 蒸腾 。 另 外 由 于 地 膜 可 以 增加 膜 内 温度 ， 进 一 
步 促 进 了 作物 根部 对 水 分 的 吸收 号。 为 了 研究 气象 因 
子 对 覆 膜 农田 蒸 散 发 的 影响 ， 研 究 了 气象 因子 与 蒸 散 


发 的 相关 性 , 如 表 3 所 示 。 气 象 因子 与 覆 膜 农田 蒸 散 
发 的 相关 性 顺序 为 净 辐 射 > 风 速 > 气温 > 相对 湿度 ， 其 
中 净 辐 射 与 蒸 散发 相关 性 最 好 ,决定 系数 R 为 0.991， 
其 次 是 风速 和 空气 温度 , ROPA 0.814 和 0.789, + 
对 湿度 与 蒸 散 发 的 相关 性 较 差 。 高 云 飞 等 所 研究 了 气 
象 因 子 对 黑河 上 游 草地 蒸 散发 的 影响 ,结果 表 明 影 响 
蒸 散发 模拟 效果 的 主要 因素 为 太阳 辐射 、 相 对 湿度 和 
Smo KAE OEO A (Phragmites australis 
Trin.) zi 38 2& BK AIA 2 rn dr AERES A 
射 、 气 温和 相对 湿度 有 较 好 的 相关 性 。 与 上 述 研 究 
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Table 3 Correlation analysis between evapotranspiration (ET) and meteorological factors of plastic film-mulched maize field 


mg 气温 净 辐 射 相对 湿度 风速 土壤 含水 量 
Item Temperature Net solar radiation Relative humidity Wind speed Soil moisture content 
R 0.789" 0.991" -0.648" 0.814" 0.093 
n 144 144 144 144 144 


** 表 示 显 著 水 平 P<0.01。** means significant correlation at 0.01 level. 


结果 相同 的 是 太阳 辐射 和 气温 与 菊 散 发 相关 性 都 良 
好 ,这 说 明太 阳 辐 射 和 气温 与 蒸 散发 的 相关 性 与 地 


膜 影响 不 大 。 然 而 上 述 研究 中 相对 湿度 与 蒸 散发 的 
相关 性 较 好 ， 而 本 试验 中 地 膜 覆 盖 下 相对 湿度 与 蒸 
散发 相关 性 较 差 ,说明 相 对 湿度 与 蒸 散 发 的 相关 性 
受 地 膜 的 影响 。 这 可 能 与 地 膜 内 相对 湿度 保持 在 较 
稳定 的 状态 有 关 , 土壤 合 水 量 与 覆 膜 农田 蒸 散 发 相 
天 性 较 差 也 可 能 与 地 膜 覆盖 有 关 。 所 以 影响 覆 膜 农 
田 蒸 散发 的 最 主要 气象 因子 是 太阳 辐射 和 气温 。 
3.2 ”影响 覆 膜 农田 土壤 蒸发 占 蒸 散 发 比例 (CE/ET) 

变化 的 因素 

土壤 蒸发 占 蒸 散发 比例 (E/ET) 的 值 反 映 了 农田 
水 分 的 利用 系数 ， 此 值 越 小 则 表示 越 多 的 水 分 通过 
植物 气孔 蒸腾 到 空气 中 ， 得 到 了 有 效 利 用 ， 在 实际 
生产 中 总 是 希望 E/ET 越 小 越 好 。 影 响 E/ET 的 因素 
除 地 膜 外 还 有 玉米 生长 环境 。 已 有 研究 表明 EET 的 
变化 主要 受 叶 面积 指数 、 气 孔 导 度 和 土壤 合 水 量 的 
影响 中 YI。 还 有 研究 发 现在 不 同时 间 尺 度 影响 
E/ET 变化 的 因素 不 尽 相同 。 在 日 变化 中 因为 叶 面积 
指数 和 土壤 合 水 量 基本 保持 不 变 ， 太阳 辐射 变化 引 
起 的 气孔 导 度 改变 则 对 E/ET 起 主要 的 控制 作用 ; 在 
生长 季 尺 度 上 ,E/ET 的 变化 受 叶 面积 指数 、 气 孔 导 
度 和 土壤 合 水 量 的 联合 影响 中 “1。 根 据 本 试验 的 结 
R, 在 生长 季 尺 度 上 叶 面 积 指数 、 太 阳 辐 射 和 饱和 
水 汽 压 改 变 引起 的 冠 层 吸收 能 量 和 气孔 导 度 的 变化 
对 E/ET 影响 较 大 ,而 在 有 干旱 胁迫 时 土壤 合 水 量 又 
成 为 限制 EET 变化 的 因素 ,这 与 前 人 的 结果 一 致 。 
通过 分 析 发 现 叶 面积 指数 主要 通过 以 下 途径 影响 
E/ET: 1) 增 加 冠 层 吸 收 的 能 量 , 减少 裸 土 吸收 的 能 量 ; 
2) 增 加 叶片 气孔 导 度 ; 3) 增 加 从 裸 土地 面 到 冠 层 的 空 
气动 力学 阻力 。 这 三 者 的 联合 使 得 叶 面积 指数 成 为 
影响 EET 在 生长 季 变 化 的 最 主要 因素 。 因 此 ， 在 日 
尺度 上 对 E/ET 影响 最 大 的 是 气孔 导 度 ， 在 生长 季 尺 
度 上 对 E/ET 影响 最 大 的 是 叶 面 积 指数 。 


4 结论 


本 试验 利用 兰州 大 学 农业 生态 系统 试验 站 2014 
年 观测 数据 ， 应 用 修订 的 S-W 双 源 模型 对 半 干 旱 区 


玉米 农田 蒸 散 量 进行 了 研究 。 以 涡 度 相关 仪 测定 的 
玉米 农田 节 散 量 为 实测 值 对 比 验证 了 S-W 模 型 模拟 
的 蒸 散 结 果 ， 分 析 了 E/ET 的 变化 太 影 响 因素 ， 并 对 
参数 进行 敏感 性 分 析 ， 得 到 以 下 结论 : 

1)S-W 模型 对 农田 蒸 散 发 的 模拟 效果 在 日 蒸 散 
量 大 于 2 mm.d- 的 晴天 和 时 畏 时 云天 气 较 好 ， 阴 雨 
天 气 模拟 结果 较 差 ， 模 型 模拟 结果 较 涡 度 相 关 的 实 
测 值 偏 高 。 影 响 履 膜 农 田 蒸 散发 的 主要 气象 因素 是 
净 辐 射 和 气温 。 

2) 履 膜 农田 土壤 蒸发 占 蒸 散发 比例 (E/ET) 在 日 
时 间 尺 度 的 变化 为 单 峰 曲线 ,在 生长 季 内 随 着 玉米 
的 生长 呈现 出 先 减 小 后 保持 稳定 的 变化 。E/ET 在 日 
时 间 尺 度 的 变化 主要 受气 孔 导 度 影响 ,在 生长 季 尺 
度 主要 受 叶 面积 指数 和 土壤 含水 量 的 共同 影响 。 

3) 敏 感性 分 析 表 明 蒸 散发 对 不 同 参数 的 敏感 性 
不 同 。 蒸 散发 及 其 组 分 对 n rn 均 较 敏感 ， 对 rn 
敏感 性 适中 ， 对 rn nj 不 敏感 。 在 应 用 S-W 模型 
模拟 农田 蒸 散 时 应 对 各 参数 合理 确定 。 
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